How close can an Inhomogeneous Universe mimic the Concordance Model? 
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Recently, spatially inhomogeneous cosmological models have been proposed as an alternative to 
the ACDM model, with the aim of reproducing the late time dynamics of the Universe without 
introducing a cosmological constant or dark energy. This paper investigates the possibility of distin- 
guishing such models from the standard ACDM using background or large scale structure data. It 
also illustrates and emphasizes the necessity of testing the Copernican principle in order to confront 
the tests of general relativity with the large scale structure. 

PACS numbers; 98.80.Cq 



I. INTRODUCTION 

The analysis of all existing astrophysical data sets, in- 
cluding the observation of the cosmic microwave back- 
ground [l|, type la supernovae (SNIa) luminosity dis- 
tance measurements Q and the large scale structure 0- 
0]) leads to the conclusion that our cosmological model 
requires the introduction of a cosmological constant, rep- 
resenting 72% of the total matter content of the Universe 
today, together with a Cold Dark Matter component 
which significantly dominates ordinary baryonic matter. 
Although such a model remains the best fit to all avail- 
able cosmological data, understanding the physical ori- 
gin of the late time acceleration of the cosmic expansion 
motivates many investigations on the basis of our cosmo- 
logical model p, 01 ■ 

As long as the Copernican principle is assumed to hold, 
there are two possibilities which could explain these ob- 
servations. Either one introduces a dark energy compo- 
nent (the cosmological constant being the simplest pos- 
sibility) or one allows for a modification of General Rela- 
tivity on astrophysical large scales (see e.g. Ref. Q for a 
recent review and Ref. Q for recent reviews on the large 
variety of dark energy models). 

The standard concordance model of the universe lies 
heavily on the assumption that the universe is spatially 
isotropic and homogeneous on large scales and can thus 
be described by a Friedmann-Lemaitre (FL) spacetime. 
Since this principle stands at the heart of our interpre- 
tation of all cosmological observations, it is important 
to ensure that it holds with a sufficient accuracy on the 
scales on which we perform our observations. Recently, 
various tests of the Copernican principle have been pro- 
posed [IMl- It was also shown that the need for dark 



energy can be suppressed if our universe is inhomoge- 
neous on scales of some fraction of the Hubble radius [Ij- 
EijilE], the underlying reason being that we actually have 
access to data only on our past null cone so that it is pos- 
sible to design spacetime geometries that will enjoy the 
same (low redshift) past-light cone structure as the one 
of a FL model without introducing a cosmological con- 
stant (T^ - [l7j (see however Ref. for some limitations) . 

In any FL model, the complete dynamics of this uni- 
verse is encoded in a single function of time, the scale 
factor a{t). This implies that all the background ob- 
servations, e.g., luminosity distance- redshift relation, an- 
gular distance, look-back time, etc., are functionals of 
the Hubble rate H{z). More importantly, in a ACDM 
model, the late time growth of density perturbation on 
sub-Hubble scales is also a function of H{z), provided we 
restrict ourselves to linear evolution. This implies that 
there exist rigidities between different set of independent 
observables ( from the background and at the level of per- 
turbation theory) that can be used to test the underlying 
hypothesis of the model 



The goal of this article is to investigate how two models 
that reproduce the same background data on the past- 
light cone can be distinguished. In Section [TTl we recall 
the main properties of the background spacetime and we 
emphasize that the time drift of the cosmological redshift 
is not only in principle but also in practice a good test of 
the Copernican principle. Section IIIII describes the per- 
turbation theory and discusses various approximations. 
In Section ITVl we describe the integration of this simpli- 
fied model in order to obtain the transfer functions on 
the past light-cone. This exercise allows us to describe 
the way one distinguishes a LTB model from a Fried- 
mann model that share the same light-cone background 
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structure. It also clearly illustrates one of the impor- 
tant limitations of most of the tests of general relativity 
on astrophysical scales since they encode the Copernican 
principle in their construction. 



II. BACKGROUND SPACETIME DYNAMICS 
A. Description of the geometry 

We consider a spatially spherically symmetric and in- 
homogeneous universe described by a Lemaitre-Tolman- 
Bondi (LTB) the metric 



1 + 2E{r) 



dr'^ + R^{r,t)dn'^ 



(2.1) 



with X = R' , using the convention R' ~ drR and R 
dtR- It is convenient to define 



2E{r) EE -fc(r)r2, 
as well as the two Hubble expansion rates 



(2.2) 



(2.3) 



The field equations for such a spacetime are given by 
M{r) 2E{r) „ „ . M' 



R^ 
i?2 



SnGpir, t) 



(2.4) 



from which it can be checked that the continuity equation 

p+(2i/i + i/||)p = (2.5) 

is satisfied. These equations can be solved parametrically 
as 



i?(t,.) = ^0'(,), 
2fc(r) 



t-tB{r) 



m{r) 
2fc3/2(r) 



(t>iv) : 



(2.6) 

where we have defined 4'iv) ~ iv ^ sinr/, rp /(S, sinhr/ — rf) 
and k = {k,r^'^ , —k) respectively for k positive, null and 
negative and defined m(r) = M{r)/r^. This solution 
involves three arbitrary functions M{r), tsir) and k{r) 
(or equivalently E{r)) but only two of them are actually 
needed since one can fix one of them by a proper choice 
of the radial coordinate r. 

To finish this description of the background spacetime, 
let us compute the expressions of i?, R' and R' that we 
will need later on. It is obvious from Eq. (|2.4I) that 



R 



M{r) 
R 



2E{r). 



(2.7) 



Then, from Eq. 



, /M' E'\ 



we deduce that 



iE' M'\ , A ■ 
(2.8) 



and then that 

IE' 



R' 



2 E 
+ t's 
1 

2R 



(2.9) 



M' 
~R 



3 E' M' 
2'e ^ M 
MR' 



M 
2R^' 



i?2 



2E' 



(2.10) 



The case of a spatially homogeneous universe is re- 
covered in the case where tsir) = 0, m =constant 
and k — constant. It follows that R{r,t) = a{t)r and 
X{r,t) — a{t) where a is the scale factor. The spacetime 
metric thus takes a FL form 



ds' = -dt^ + a\t) [dx' + fK{x)dn^] 



(2.11) 



with fxix) ^ (shi-X, X, sinh x) depending on the sign of 
K. 



B. Light cone equation 

Most of our analysis focuses on what is actually ob- 
served on our past light-cone, which can be defined by 
solving the null geodesic equation. Given a null geodesic 
with tangent vector fc^, the redshift of any object is de- 
fined by 



re 



(2.12) 



where is the tangent vector to the matter worldlines. 
It follows that the geodesic equation k^^V ^k'^ = for a 
null- vector reduces to 



dt 
d^ 

dr 



1 



{l + z)H{ 
v/1 + 2i?(r) 



(2.13) 
(2.14) 



dz (1 + z)R' 

This defines our past lightcone that we shall denote by 



{r = r,{z),t^ U{z), n(0) = 0,t,(0) = h}. 



The redshift is the observational radial coordinate which 
can be expressed as a look-back time or distance, assum- 
ing a cosmological model [2l|. One can indeed use ei- 
ther r, t on z as integration variable, the important point 
being that only 2 of the three quantities (t*,r*,z) are 
independent. 

When evaluated along the past light-cone, the function 
R{r, t) is related to the angular distance Da by 



7^o(z) = R[U{z),r,{z)] ^ Da{z) . 



(2.15) 



This observational relation can be used to set one con- 
straints on the two arbitrary functions of the LTB space- 
time. It is useful to redefine the functions introduced 
previously when evaluated on the past-light cone as 
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R[U,r4, nQi = R'[t,,r,], 



(2.16) 
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and 



We thus impose that 



7^ll = R'[U,r^ 



(2.17) 



to be considered either as functions of z or 

We choose the radial distance along the hght cone 
and fix it by imposing 



7^n 



^/l + 2E 



1, 



(2.18) 



which simphfies the past light-cone equation reducing it 
to 

^ ^ (l + z)7eii(rH.) 
dr* ' dr, ^1 + 2S(n) ■ 



C. Reconstruction procedure of a LTB geometry 

1. General procedure 



With the choice of coordinates ()2.18p the derivative of 
TZ along the past light-cone is given by 



d7^ 



= 7^ol-7^ 
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(2.20) 



and Eqs. (|2.8ll2.9p can be rewritten as 



M' 

[7^o -{t- tB)nw] — 



-(^-^B)7^lo 



'e 



'TZint', 



M M' 



2nl M 



7^ol = ^/l + 2E, (2.21) 
El. 

'e 



3, , M 
7;{t~tB)- 



2'-" ^"'2nl 



M 



2TZ, 



t'n — Till 



(2.22) 



Hence, provided Da{z) is known from observation, we 
fix 7?,o and then Eqs. (|2. 19112. 22| give 4 equations for the 
6 unknown functions {t,z,M,E,tB,'R-) of r*. One thus 
need to fix 2 conditions to completely specify the model. 

This approach was investigated in Ref. [la| who de- 
signed various LTB solution sharing the same Da{z) re- 
lation as a FL spacetime by imposing either = or 
k = constant. We follow the same line and consider LTB 
models that reproduce the angular distance-redshift re- 
lation of a fiducial standard fiat ACDM model. For such 
a FL universe 



^0 



= f7m0(l+z)3 + f7A0, 



(2.23) 



with fl\Q = 1 — rimo- Hq is the Hubble parameter eval- 
uated today and the density parameters are defined by 
f^mo — STrGpmo/S-ffo ™d flAo = A/3Hq. The angular 
diameter distance is then given by 



Da{z) = 



1 



du 



Ho{l + z)Jo ^r!„,o(l + u)3 + n 



(2.24) 



n{z) = DA{z), 



(2.25) 



so that both cosmological models enjoy the same angular 
distance-redshift relation (and thus luminosity distance- 
redshift relation, as can be observationally determined 
from SNIa observations). 



2. A model mimicking the background light-cone dynamics 
of a FL universe 

To completely specify the model, we need to impose 
another condition on the LTB model. In order to con- 
struct a model that would mimic a fiducial fiat A-CDM as 
close as possible, we follow Ref. Q and further assume 
that the LTB matter energy density distribution along 
the light-cone, p{z), matches the observed mass density 
as a function of redshift as determined in the fiducial 
ACDM model, that is 



8ttGp{z) = 8ttGpfl{z) = 3n,^oH^{l + 



(2.26) 



This assumption allows to determine r*(z) (T^. Il7|. and 
in the particular case of a fiat ACDM, this implies [l^l 



dr* 
d^ 



{1 + z)Hfl{z)' 



(2.27) 



AO 



Comparing with Eq. (|2.14p . this means that, under the 
choice (|2.18p for the radial coordinate, H^^{z) = Hp^^z). 

We then follow the reconstruction procedure described 
in Refs. [H, [T^l and our result agrees with those pub- 
lished in these works. Figure [T] depicts the function m(r) 
and k{r) obtained from this reconstruction. 

Eqs. (|2.7II2.9P imply that when r ^ 0, R/r 1 while 
TO — ^ finio-ffo ^i^d k — >■ (rimo — ^)Hq. The last free func- 
tion i_B(r) can be reconstructed from Eq. (|2.9p but we 
will not need it explicitly in the following. Imposing that 
is(0) = 0, we can compute to from (|2.6p . after evaluating 



D. Distinguishing the two models 

The LTB model defined in the previous section was 
designed to strictly mimic the luminosity (or angular) 
distance redshift relation and the mass density-redshift 
relation of the fiducial ACDM model. 

In order to distinguish these two models, one needs to 
find another independent observable quantity. It was re- 
cently pointed out that one can extract some information 
of the dynamics off the light-cone by considering the time 
drift of the cosmological redshifts 'lOi] . While advocated 
as a test of the Copernican principle its amplitude in a 
non-FL model was not estimated. This can easily be ob- 
tained for the LTB model under investigation. We first 



FIG. 1: Reconstruction of the function m(r) (black, solid line) 
and k{r) (blue dashed line) entering the definition of the LTB 
geometry for a spacetime reproducing both Da (2) and pmiz) 
on the past light-cone. The light dotted lines correspond to 
m(0) and fc(0). 



remind that the time drift of cosmological redshift in a 
LTB universe takes the form [10, [2311 



(2.28) 



which generalized the original FL-expression [2J| . Indeed 
in a FL-model H± ~ H\\ ^ H so that i derives from 
H(z), as any other background observations. However, 
in our LTB-model H± ^ = Hp^. This was used to 
demonstrate that z[z] allows to fully close the reconstruc- 
tion system without resorting to making assumption on 
the matter energy density profile along the light cone. 

Figure [5] compares the expected time drifts of the cos- 
mological redshift for the LTB- and FL-models. Az — 
zAtohs has a typical amplitude of order 10^^ on a time 
scale of Aiobs = 20 yr, for a source at redshift z ^ A. 
This measurement is impossible with pre sent-day facili- 
ties. However, it was recently revisited (25| in the context 
of ELT, arguing they could measure velocity shifts of or- 
der 110 cm/s over a 10 yr period from the observation 
of the Lyman-a forest. It is one of the science drivers 
in design of the CODEX ultrastable spectrograph [2^ 
for the future European ELT. Indeed, many effects, such 
as proper motion of the sources, local gravitational po- 
tential, or acceleration of the Sun may contribute to the 
time drift of the redshift. It was shown [13, HI]) however, 
that these contributions can be brought to a 0.1% level so 
that the cosmological redshift is actually measured. The 
data points and error bars of Fig. [2] follows the forecast 
of Ref. ^2d\ . Our analysis confirms the recent analysis by 
Ref. [29] , which also suggests to use the cosmic parallax. 



E. Discussion 

The result of Fig. [5] demonstrates that the information 
off the light-cone allows to distinguish a A-CDM from a 



FIG. 2; Time drift of the cosmological redshift for the stan- 
dard A-CDM model (black, solid line) and a LTB-model (blue, 
dashed line) designed to share the same observational relation 
on the past light-cone. The data points follow the estimates 
of Ref. for a CODEX-like spectrograph on an ELT. 



LTB-model specially designed to have the same luminos- 
ity (or angular) distance redshift relation and the same 
mass density-redshift relation as the a A-CDM. In the 
case of a dark energy model both Dl{z) and z are mod- 
ified, and little insight is gained on the equation of state 
from adding the new information on z 3^. The case 
of the large scale inhomogeneity turns out to be differ- 
ent and the i observation would provide, when available, 
an interesting extra-test of these models that cannot be 
performed otherwise. 

Indeed, Fig. [2] was obtained by forcing Da{z) to match 
the ACDM prediction up to 2: ~ 4 while SNIa data ex- 
tend roughly to z ^ 1.6. Following e.g. Ref. [Slj, one can 
try to design density profiles such that the LTB model re- 
produces the ACDM-D^(2;) at low redshift and becomes 
homogeneous on large scales. Note that reproducing p 
and Da of a flat- ACDM imposes that H±{z) = Hfl{z). 
Now, assuming that it fits zpL and p at high redshift im- 
plies that H± = — HpL- While attractive, such mod- 
els seem however difficult to construct since they usually 
require that M' < between these two regimes. 

This example demonstrates the complementarity of 
these two observables since they concern two domains 
of redshift. To go further in distinguishing such a model 
from its FL-twin, we shall now consider the influence of 
the large scale inhomogeneity on the growth of large scale 
structures. 



III. EVOLUTION OF DENSITY 
PERTURBATION IN A LTB UNIVERSE 

A general study of the perturbation theory around 
a LTB background was performed using a coordinate 
based approach in Ref. [33], and the general features of 
the grow th of density perturbations were discussed in 
Ref. [3^. The goal of this section is to investigate the 
evolution of the density contrast using the 1-1-1-1-2 for- 
malism, and to obtain an approximation for the evolu- 
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tion equations. We also discuss how the density contrast 
variable introduced here can be related to observations. 



A. General formalism 

1. 1+3 formalism 

In the 1+3 covariant approach [34], in which one intro- 
duces the worldline tangent vector {u°'Ua ~ —1), one 
first introduces the projection tensor = S"-^ + u°'ui,. 
This projection tensor defines two derivatives for any ten- 
sor T°'hi the covariant time derivative along the funda- 
mental worldlines 



(3.1) 



and the fully projected covariant derivative I?e via 

D,T\ = h''fh%h\VrTfg , (3.2) 

where we fully project on all free indices. 

The projection tensor hat allows for any 4- vector to be 
split into a scalar part parallel to m° and a 3- vector part 
orthogonal to u". Similarly, any second rank tensor may 
be covariantly and irreducibly split into scalar, vector and 
projected, symmetric, trace-free (PSTF) 3-tensor parts. 

We can also split the covariant derivative of Ua into its 
irreducible parts as 

VaUf, = - Uailb + DaUb (3.3) 
= -UaAb + lQhab + (Jab + ^ab ■ (3.4) 

This uniquely defines the following kinematic quanti- 
ties: Ab = Ub is the acceleration, the trace Q = DaU°- 
is the (volume) rate of expansion of the fiuid (with 
H =^ Q/'S the Hubble parameter), aab = D(^aUb) is the 
trace-free symmetric rate of shear tensor describing the 
rate of distortion of the matter flow, and uJab = £'[0^6] is 
the skew-symmetric vorticity tensor, describing the rota- 
tion of the matter relative to a non-rotating frame. It is 
also useful define a scale factor a along the fundamental 
worldlines via 



1 



e. 



(3.5) 



2. 1+1+2 formalism 

We now employ the 1-f 1-1-2 formalism fss'], that builds 
on the 1-1-3 formalism by allowing a further "spatial" slic- 
ing with respect to a spacelike unit vector field n° , which 
is orthogonal to the timelike 4- velocity vector {u°-na — 0). 
The 1+3 projection tensor hj' combined with gives 
rise to a new projection tensor iV^^, which projects vec- 
tors and tensors onto the 2-surfaces orthogonal to 71° and 
u"' that are referred to as the 'sheets': 



Analogously to the 1+3 formalism, Na^ defines two new 
derivatives for any object T... ": 

iT\y = n''D,T\, (3.7) 
d,T\ = N'^fN\N\DrTfa. (3.8) 

The '-derivative is the (spatial) derivative along the vec- 
tor field n° in the surfaces orthogonal to analogous to 
the time derivative defined in (j3.ip , and the c?e-derivative 
is the projected derivative on the sheet analogous to the 
projected derivative De defined in (|3.2p ^. See Ref. [s^ 
for a detailed presentation of the 1+1+2 formalism and 
the appendix for a summary of the key equations. 



3. The LTB case 

For a LTB spacetime all vectors and tensors as well 
as vorticity, acceleration and the magnetic part of the 
Weyl tensor vanish at background level. It can be fur- 
ther shown that the expansion, the only non-vanishing 
component of the shear E = n^n^a^i, and of the electric 
part of the Weyl tensor 8 — n^n'^E^i, are given by 



e = 2i7 1 



and 



E = 



SnGp M 



(3.9) 



(3.10) 



We depict their evolution on the past light-cone in Fig.|3l 
and they satisfy the following general evolution equa- 
tions: 



e = 
± = 

£ = 



-ie^-4.Gp-^I]^ 



1 



e - -E ) f - inGpT,. 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



B. The density perturbation equations 

We can decompose the standard dimensionless normal- 
ized density gradient Va = a hab^''p/p into a part Aa 
that lies in the sheet and a part perpendicular to the 
sheet: 

T. P' daP _ p' , - 

Va^a—Ua + a = a— na + Ao , (3.15) 

P PR 



N"-b = h'^b - n^fib = g'^b + u'^Ub - n'^Ub . 



(3.6) 



^ Note that the '-derivative was denoted by a hat Q, and the de 
was denoted Se in Ref. 13511 . 
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to 
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FIG. 3: Evolution of the expansion O (black, solid line), the 
shear E (blue dashed line) and the electric part of the Weyl 
tensor £ (red dotted line) on the past light-cone for the LTB 
model defined in Section Hi CI 



where we recall that p' = n^Dhp and daP = N^aDtP- 

It is clear that is gauge-invariant as it vanishes in 
the background LTB model. To see more clearly what 
this variable physically represents, let us define an in- 
finitesimal vector 6xa = N''a Sxt connecting two neigh- 
boring points on the sheet. It follows that the density 
evaluated at two neighboring points are related by 



p{Xa + SXa) = piXa) + daP Sx"" . 



(3.16) 



In this way we can define a gauge-invariant density con- 
trast S: 

5 ^ Pi^<^+S-^) -PM ^ l^J,-a . (3.17) 

p(Xa) a 
Defining the auxiliary perturbation variables in the sheet 

Se = -eaSx"", = -TaSx"", 5£ = -SaSx"" 
a a a 



where 



9a = a daQ , Ta = a doT, , Sa = a da£ , (3.19) 



we can derive a set of 4 gauge invariant perturbation 
equations (neglecting the magnetic part of the Weyl ten- 
sor) 



(3.20) 



(3.21) 



S = -SQ + !-CaSx'' 
P 

[6QY = -^e(5e-47rGp<5-3S(5E 

-e'c-aSx"" 

(5E)- = ^-^e - (51] - ^E(5e - 5E 
-Y.'C-aSx'' 

-£5<d - AnGpYsS - E'C-aSx'' , (3.23) 



(3.22) 



where Cg = T,a + ota- The equation for the connecting 
vector 5x°' is only needed to background order. 



2 V3 



(3.24) 



The set of equations 



provides a complete de- 
scription of the evolution of the gauge invariant density 
contrast in situations where the magnetic part of the 
Weyl tensor can be neglected. This 'silent" approxima- 
tion can be thought of as neglecting the coupling be- 
tween density perturbations and gravitational waves in 
pressure-free models and has been widely studied in the 
context of non-linear gravitational collapse [53 |. 

The above equations can be solved iteratively as fol- 
lows. Expanding the perturbation variables in the form 
= 5Xo -f pjx'', (with i = 1 ... 4 and the short-hand 
notation 5X^ = 6 etc.), we obtain a set of homogeneous 
equations for the leading part of the perturbations 



^0 



-590 



(5e)b = --QSQo-^7TGpS~3ESEo 



(3.25) 
(3.26) 



(,5E)b = (-^e-E^5Eo-^E5eo-<5fo (3.27) 



-£S& - 47rGpE(5o ■, 



(3.28) 



with the corrections evolving according to equations of 
the form 

/a^-^Qe-E^ji^ + a^Ca (3.29) 



with a* = (p'/p, — 0', — E', — f ). The homogeneous 
equations can then be combined to give a pair of cou- 
pled second order equations 



So + -e(5o-47rGp(5o = 3E(5Eo, 
(<5E)o + (^e-iE)(5E)o 



(3.30) 
(3.31) 



20 2 4 1 

-7:Gp + -GE - -e^ --£ + 5E2 ) JEq 



4,2 2 \ ■ 20 

-E2 + -f + -Eej5o + y^Gp,5o. 



The first one is similar to the standard one in FL models, 
but with a source term proportional to the background 
and perturbed shear so that it is coupled to a second 
differential equation for ^Eq. 

In the cases where the above approximations are not 
valid, one needs to consider the considerably more com- 
plex full set of l+H-2 equations that is described in the 
appendix. 
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IV. TRANSFER FUNCTIONS AND INITIAL 
CONDITIONS 

A. Back to Friedmann-Lemaitre spacetime 

Before investigating the growth of structure in a LTB 
universe, let us reeah the relations between the back- 
ground dynamics and the growth of large scale structure 
in a fiat ACDM model. 

Since the cosmological constant develops no perturba- 
tion, and since the growth of the density perturbations 
of the pressureless CDM component is dictated by gen- 
eral relativity, it implies that in the linear regime, the 
evolution equations reduce to the two equations 

S = -(59, SQ = -2HpLSe - AnGpS , (4.1) 

that is to 



D + 2HD - 4:TTGpn,D = 0, 



(4.2) 



where the CDM density contrast has been decomposed 
as S{x,t) — D+{t)e+{x) -I- D_(i)e_(x), e± encoding the 
initial conditions. This equation can be recast in terms 
of a as time variable (sgI . Ist} as 



D" 



dlniJ 
da 



mO 



D. 



(4.3) 



Since D_ = _ff is a solution, the growing mode is given 

by 



mO 



(1 + z')dz' 

WWW' 



(4.4) 



In our particular case, this can be integrated analytically 

as 



D+{z) cx 2F1 



1 11 



1 - 



mO 



mO 



mO 



1/3 



1 



1 



[l + zf 



1 



(4.5) 



This implies that if H{z) is known from background 
observations, such as SNIa, then -D+(z) is fixed and is 
thus not an independent quantity in this framework. 
There is a rigidity between the expansion history of the 
background and the growth rate 3^ [s^ . This has at- 
tracted attention since this offers a test of the ACDM 
model m, SO]. Moreover, Eq. dUI]) implies that = 
SQ/H is related to the density contrast by 

d = -f{a)6 (4.6) 

where the function / can be parameterized as [ill, l42l | 

d\nD. 



din a 



(4.7) 



for the growing mode. Then, if general relativity is not 
modified, the index 7 can be computed once H{z), or 



equivalently the dark energy equation of state, is known 
and it was shown [4^ that 7 = 0.55 in the case of 
a ACDM. Since 7 can be measured from galaxy redshift 
surveys (iG^ , thanks to the redshift distortion [4^ , it was 
argued that the value of the parameter 7 offers a test 
of general relativity [39[ . 



B. Structure of the LTB growth rate equations 

In a LTB spacetime, and under the approximations 
discussed in Section IIIIl we need to solve a set of 4 dif- 
ferential equations that can be recast as 



(4.8) 



where Xi — {S, 6Q, JS, 6£) and where the matrix Mij de- 
pends only on r and t through the background quantities 
O, S, p and £. It follows that the general solution is of 
the form 



Xi{r,t,0,(p) = Tij (r, t)Xj (r. Unit ,0,Lp) 



(4.9) 



where the angular dependence stems only from the ini- 
tial conditions so that the transfer functions depends 
only on r and t. When interested in observations such 
as weak-lensing, one integrates along the light-cone so 
that only (z) = Tij [r* (z), t,(z)] is actually needed (see 
Ref. [33|]). 



1. Integrating the perturbation equations 

In order to compute the transfer functions, we proceed 
as follows (see Fig. S]). 

• Using the reconstruction along the past light-cone, 
wc integrate for each z the background equations 
(|3. 11113.1^ from t = to ~~ r*(z) toward the interior 
of the past light-cone at r = r*(z) constant. This 
provides the background quantities p{t,r), Q(t,r), 
E(t, r), and £(t, r) for this particular r in agreement 
with the constraint that R — Da and p = ppL on 
the past light-cone. 

• We then solve the perturbation equations with this 
background functions from an initial time iinit up to 
i*(z) considering the four sets of initial conditions: 



init ; 



a = 1 . . .4. 



(4.10) 



For each set a, we obtain 4 transfer functions Ti^. 

Such an integration procedure is made possible because 
the background equations involve no gradient term, so 
that each shell of constant r evolves independently and 
because, in the silent approximation used in this work, 
the perturbation equations enjoy the same property. In 
this approximation, the spatial inhomogeneity of the 
background spacetime reflects itself only on the fact that 
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a structure observed at a redshift z had a growth history 
along the sheh r = r»(z) which is difTerent from the other 
sheh. This is a major difference with the FL situation in 
which, as soon as we are deahng with a pressureless fluid, 
all the structures have the same growth rate from z to 0, 
independently of their growth rate at higher redshift. 




FIG. 4: Summary of the integration procedure. All quanti- 
ties are defined in the text. From the light-cone background 
quantities one integrates the background equations toward 
the interior of the light-cone at constant r (1). Then, once 
an initial hypersurface has been chosen, one integrates the 
perturbation equations up to the light-cone (2). 

First, this procedure can be tested on the FL model 
for which E = so that ((5,(50) decouples from (5E,(5£). 
Focusing on the matter density contrast S, we compute 
the two transfer functions Tpp and Tpg corresponding re- 
spectively to the initial conditions ((5,(50) = (1,0) and 
((5,(50) — (0,1). None of these transfer functions cor- 
responds to the growing mode (|4.4p that is obtained 
from the initial conditions ((5,(50) = (1, — /+(<init)-ffinit) 
so that 

D+{z) = Tpp{z) ~ f+{Unit)HinitTpe{z). (4.11) 

Figure [5] shows that this is actually verified numeri- 
cally. Note that even though we reduce the dimen- 
sion of the space of initial conditions by picking up 
the growing mode, we still need the two transfer func- 
tions. Indeed one can also check that S9+{z) = Tgp{z) — 

2. Initial conditions 

To compute the transfer functions of the LTB model, 
we need to specify an initial spatial hypersurface. Any 
3-dimensional spacelike hypersurface {r = tinit('')} is a 




- 

^ " ~ 7 : : : r : : : r : 

0.2 0.4 0.6 0.8 1.0 

1/(1 +Z) 

FIG. 5: The two transfer functions Tpp (dashed blue line) 
and Tpe (dotted green line) for a ACDM. We check that Tpp — 
f+{zinit)H {zinit)Tpe (solid black line) is the growing mode D+ 
(dots). 



priori possible but three choices can be argued to be nat- 
ural: (i) a constant time hypersurface, (ii) a constant 
density hypersurface or (iii) a constant t — tB(r), i.e. a 
constant proper time after the big bang. Indeed, in a 
FL model, these three possibilities reduce to the same 
hypersurface. 

For the purpose of the illustration, we decide to set 
the initial conditions on a constant time hypersurface, 
as in Rcf. 33J, but we cannot justify this choice further 
here. Then, applying the procedure described in Sec- 
tion IIVB 1[ we obtain the transfer functions. Figure [S] 
describes the four transfer functions needed to compute 
the density contrast on the past light-cone. 




12 3 4 

z 



FIG. 6: The four transfer functions Tpp (sohd black line) and 
TpB (dashed blue line), Tpa (dot dashed red line) and Tpg 
(dotted green line) for the LTB model. 



3. Comparison of the two models 

With the previous prescription, the two transfer func- 
tions Tpp and Tpg for the two models look similar in shape 
and amplitude. 
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To compare the growth rate, one needs to know the 
initial conditions, or at least the relative initial value of 
5, 5Q, STj and 6£. As we have seen in Section flVB H a 
particular choice of {S, SQ) allows to pick up the growing 
mode. We have no indication of the linear combination 
of the transfer functions that are related to the growing 
mode in the LTB case. In order to compare the growth 
rate, we assume that at early time the universe was well 
described by a FL-model and that the void evolves at 
lower redshift. In such a case, the initial conditions can 
be set by their FL analog for the growing mode. From 
Eqs. (|A28|) and (|A31|) . we deduce that 

-9(5Einit = -(87rG(5pinit), -^SSinit = -^-^Spinit 

(4.12) 

and for the growing mode 

<5ei„it = -/+(init)il(init)5i„it. (4.13) 
That would imply that S = DY^^Si with 



D 



LTB 



(z) = Tpp{z) - f+{mit)H{mit)Tpeiz) (4.14) 



4 2 
——Tps{z) - -Tpsiz) 



It corresponds to the growing mode of the density pertur- 
bation if the universe has evolved from a FL-phase and if 
the density perturbations had time to reach the FL grow- 
ing mode before the effect of the void on the evolution 
of the perturbations starts being non-negligible. Indeed, 
we have no proof that it is the growing mode of the LTB- 
system. Note that the 4 transfer functions are needed to 
described the evolution of the density, even though the 
initial conditions can be reduced to the single random 
variable (5i„it(r, tinit, 0, ip). 

Fig. [7] compares this solution to the growing mode of 
various FL-models. All models are normalized at high 
redshift, and we can see that the LTB-model has more 
structure at small redshift and that its growth rate is 
qualitatively similar to the one of a closed FL-model. 
We can also note that it seems impossible to find a lin- 
ear combination of the 4 transfer functions that would 
mimic the growth rate of the ACDM model. This is an 
interesting conclusion, given our ignorance of the initial 
conditions. Indeed, in general, we expect the initial con- 
ditions to be not simply a linear combination but of the 
form 

<56mit = f(d{r)5init{r, 0, (fi), 5Sinit = h{r)Sinit{r, 0, (fi), 

(4.15) 

and 



(4.16) 



that is involving three arbitrary functions of r. Unless 
we have theoretical constraints of /e, /e and /f , the 
growth rate (j4.14p can be tuned at will. It follows that 
the determination of the growth rate is more likely to 



teach us on the initial conditions than to be used as an 
extra-data for the reconstruction program. 

As emphasized, in the ACDM, the growth rate and 
the background dynamics are not independent so that 
we have relation such as Eq. (|4.7I) . For the ACDM, 7 = 
0.55 and slight deviation from this value are expected 
for dark energy models in which general relativity is not 
modified [391 • The index 



e(a) = [fJ,„(a)]" 



dlnD 
din a 



1, 



(4.17) 



introduced in Ref. 40] should not deviate significantly 
from 0. 7 and Qrn are determined from background ob- 
servations and thus coincide with their ACDM values. It 
was shown that e can typically vary between 0.05 and 
0.25 for modification of general relativity of the f{R)- 
class poj . Using the numerical solution corresponding to 
Fig. [71 we estimate that e can reach 0.1. The deviation 
arises from the fact that the perturbation equations in- 
volve the shear and the electric part of the Weyl tensor 
and indeed not from a deviation from general relativity. 
This illustrates that it is indeed important to ensure the 
validity of the Copernican principle when applying tests 
of general relativity based on the large scale structures 
since they usually implicitly assume its validity. 



Q 

< 



FIG. 7: Comparison of tlie density growing mode for a flat 
ACBM (blue dashed line) an Einsein-de Sitter model (red 
dotted line) and our LTB model (black solid line) assuming 
the FL-like initial conditions. It can be mimicked by a closed 
FL (green dotted line). 

To finish, let us also emphasize that any observable, 
O say, entangles the properties of the transfer functions 
and of the initial conditions. In general, it can be written 
as an integral along the line of sight as 

0(z,n")= f w{z,z')6[x''{z',n'')]dz', 
Jo 

where S stands for some combination of the perturbation 
variables and w for some window function that depends 
on the observable (e.g. for weak lensing in a FL universe S 
will correspond to twice the gravitational potential and 
w can be expressed in terms of angular distances). In 
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a Friedmann-Lemaitre universe, as long as we focus on 
the evolution of dark matter, we have that S{x'^,t) ~ 
T{t)dinit{x^) so that the correlation function of O takes 
the form 

{0{z,n1)0{z,n^)) = f dz' f dz"w{z, z')w{z, z") 
Jo Jo 

T(z')T(z")6nit ( ^r2(z') + r2(z") - 2r{z')r{z") cos0) , 

with cos 9 = ni.n2 and where the correlation function 
of the initial condition, ^init is a function of \x\ — ajjl 
only because of isotropy and homogeneity. This integral 
is usually easily evaluated in Fourier space. On small 
scales, the initial conditions decouple from the evolution 
mainly because the power is mainly carried by modes 
perpendicular to the line of sight (which is at the basis 
of the Limber approximation) but the initial power spec- 
trum has to be such as this property holds, which is the 
case for an almost scale invariant power spectrum. In 
a LTB universe, there are two main differences with the 
FL case. First the transfer functions depend on r and 
t so that S{x'-,t) = T(r, t)(5init(a^*) and second the initial 
correlation dunction is a function ^init{r{z'),r{z"),cos9) 

{0{z,nl)0{z,n^)) ^ J dz' J dz"w{z, z')w{z, z") 

T[r{z'), i(z')]r[r(z"), t(z")]6nit [r(z'), K^"), cos 9]. 

It follows that the spatial dependency of the transfer 
functions mixes with the one of the initial conditions and 
that we cannot ensure a priori that a condition similar 
to the Limber property will exist. 

Indeed, once a theory of the initial conditions compat- 
ible with the origin of a large scale inhomogeneity exist, 
one can use our analysis to compute the angular power 
spectrum of the observable O. While being focused on 
the properties of the transfer functions, our work however 
show that if ^initi is identical as in a FL universe then the 
effect of the evolution can lead to significant effect. 

V. DISCUSSION 

In this article, we have used a LTB model that mim- 
ics a FL model on the past light-cone in order to shown 
that these two models can still be distinguished by back- 
ground observations that encode information "off" the 
light-cone, as for instance the time drift of cosmological 
redshift. While such an observation was advocated as 
a test of the Copernican principle in Ref. 1^ , the am- 
plitude of the difference with the ACDM prediction had 
not been estimated. We have shown that it can be sig- 
nificant and that it can allow to exclude a large class of 
LTB models even though they look similar to FL models 
at the background level. 

Then, we investigated the information that can be ex- 
tracted from large scale structure. Assuming that the 
curl of the magnetic part of the Weyl tensor can be ne- 
glected, we have shown that one can extract a closed 



system of 4 perturbation equations. Our derivation clar- 
ifies the link with the gauge invariant variables of the 
1 + 1 + 2 formalism. On small angular scales, we argued 
that one can make a silent approximation. Under such 
conditions, our set of equations for scalar perturbations 
reduce to the ones used in Ref. [ssl ]. 

We have detailed a procedure to compute the transfer 
functions of this system of perturbation equation, once 
the background data is known on the past light-cone and 
we explained how the corrections at first order in the 
connecting vector can be obtained. We emphasized the 
difficulty in determining a set of the initial conditions, 
i.e. to define the hypersurface on which they are defined 
and the relation that could exist between the different 
modes (in particular to extract the growing mode). 

Indeed, our discussion of the perturbation dynamics is 
more illustrative than quantitative for two reasons. The 
first one is mainly technical since we have used the silent 
approximation instead of solving the full set of 1+1+2 
equations. This can only be achieved numerically and is 
left for further investigation. The second is related to an 
intrinsic limitation concerning our ignorance of the ini- 
tial conditions: (1) we do not know the hypersurface on 
which they have to be specified and (2) while we have 
been able to compute the transfer function, we are not 
able to specify their exact combination that corresponds 
to the growing mode. While in a FL universe, the late be- 
haviour of the growth function of dark matter reduces to 
a function of time, or equivalently of redshift, it is in the 
LTB case a function of t and r that reduces to a function 
of z on the light-cone. This implies a further limitation 
since the large scale structure properties strongly entan- 
gle the evolution and the initial conditions. It is thus dif- 
ficult to close the reconstruction program unless we have 
some constraints on the initial condition. Nevertheless, 
our investigation shows that we cannot find a linear com- 
bination (with constant coefficients) of the transfer func- 
tions that reproduce the growing mode of the ACDM. It 
also illustrates explicitely one limitation of the tests of 
general relativity based on the large scale structure (e.g. 
the 7-index is difficult to generalize to a LTB spactime 
because Eqs. (|4. 6114.71) have to be reconsidered). 

This emphasizes the importance to test the Copernican 
principle to validate the test of general relativity based 
on the properties of large scale structures and illustrates 
the effect of the assumption on the large scale geometry 
of our universe on these tests. In that respect both the 
time drift of the cosmological redshifts and the growth of 
density perturbations give access to the spacetime struc- 
ture beyond the light-cone and are thus key observations 
for our understanding of the geometry of our universe. 
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Appendix A: Full set of perturbation equations in the 1+1+2 formalism 

If the magnetic part of the Weyl tensor or the source term (Sa + aa)(5" cannot be neglected, the linear perturbations 
around a LTB background are no longer described by equations (13.20113. 2^ . and the complete set of first order equations 
below will successively be required to obtain a closed system of linear perturbations. 

Let us first recall that in the 1+1+2 framework any 3-vector or PSTF 3-tensor Tat can be irreducibly decomposed 
with respect to into scalar, vector and tensor parts: 

Tab T {nam - ^Nab) + 2T(anb) + Tab, (Al) 

where 

r = n'^V^Tab = -N'^'Tab, 

Ta = Nju'^Ttc, 

Tab = (n^^-N^)" ~\NabN''')Td. (A2) 

The shear aab e.g. will be decomposed into the scalar part S = n°'n^aab, vector part Eq — N^rf'crbc and tensor part 

^ab - [n^^^N^^'^^ \NabN-''') 0,d. 

We may also decompose the covariant derivative of n° orthogonal to in analogy to the 1+3 decomposition (13. 3|) : 

D^rifc = TlaCfc + i^iVafc + ^ffafc + Cab, (A3) 

where 



Oq = n^DcHa = n'^, (A4) 

^ = dan'', (A5) 

e = ^e'^'dam, (A6) 

Cab = d^anb}- (A7) 



Along the direction n", (j) represents the sheet expansion, Cab is the shear of n", a" its acceleration, while ^ represents 
a 'twisting' of the sheet. For the LTB background. Cab, n"" and ^ are first order quantities. Furthermore, we restrict 
ourselves to perfect fluid perturbations where pressure, 4-acceleration and cosmological constant vanish at all orders. 



1. Source evolution 



The source evolution is then governed by the propagation equations 



Sa = - 0e + Y.a - h^Oia " £a (A8) 

E-a + le^V!^ = \eabd''U + \ebcd''U\ - \ {%i^Gp - |f ) - f + (|S - 9) - i^Eab^'' (A9) 

S{a6} = - (|e + iE) Eab - fab (AlO) 

£{ab} - £c{aK}" = -Sciad^Hb} - \ i&nGp + 3£) ^ab - (6 + fS) Eab + i0£c{aHfc}' (All) 

n = -sabd^E^ -iS^E + {id + p:)n (A12) 

Ha-iea'ffc - -leabd'E^\ebcd'E^, + lEeabo' + \<i>eabE'+{l^-Q)Ha (A13) 

■H{ab} + £c{affc}" = ec{arf'fb} + ^claCb}' - 5'^ec{afb}^ " (6 + f S) Hab (A14) 

Oa-a'a = \<i)aa + {\Q + Y)aa-\<l>'i^a + £ab'H^ (A15) 

= -_(20_2) i0 + d„a° (A16) 

1 - (iS-ie)^+i£,bd"a^ + iH, (A17) 

C{ab} = {^"^ - Cab-\(t>'^ab + d{aab} - £c{aT-L^ , (Al8) 



where we use curly brackets to denote the PSTF with respect to n° part of a tensor. The last equation may not be 
required up to first order, but just at background level, since it only seems to appear coupled to first order terms. 
Note that there is no propagation equation for Oia ■ 
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2. Auxiliary First Order equations 

Auxiliary first order propagation equations along u° are given by 



ie-E = (ie + S) +|7rG'/9 + f (A19) 

-e = ie2 + |l]2+47rGp (A20) 

E = -(|e + iE)E-f (A21) 

p = -Op (A22) 

£ = (|E-6l)£:-47^GpE + eabrf"■H^ (A23) 



and the first order propagation equations along n° are 

(/,' = -102 _ (10 ^. 2) (|e - E) - f TrGp - £ + daa" (A24) 
^' - + ^eafcrf'a" (A25) 

C{ab} = -<^Ca6 + rf{aab} + + E) Eab - £,6 (A26) 

E'-|e' ^ -f(/)E-daE'' (A27) 
E;, = idaE+|dae-|(/.E, -fEoa-d'Eah (A28) 

^{ab} = C?{aSfc} - 50Eah + §ECq& - £c{aH^ (A29) 

£"-f7rGp' = --daS" - l(t>£ (A30) 

= idaf + f TrGdaP - dXb - fla - a - EeabH*" (A31) 

n' = -daW-^m (A32) 

= idaH - d''Hab - |f eafcE^ + i&abf ' - (A33) 

Appendix B: Weak lensing for a central observer 

One of the key observation to extract information about the growth rate of the large scale structure is weak lensing, 
in particular using future tomographic survey. We recall the Sachs equation [s^, 113, El] , which is the central equation 
governing gravitational lensing and then consider it in a LTB universe. 

1. Sachs equation 

The tangent vector to a null geodesic is fc° — dx°/dA and satisfied {1 + z) = k"-Ua if we choose the value of the 
affine parameter A such that k'^Ua — 1 today. It follows that it can be decomposed as 

r = -(l + z)(M"+n"). (Bl) 

is the spatial direction of the photon and, in the particular case in which the observer is at the center of the 
spherically symmetric spacetime, the null geodesies are radial and n° reduces to the radial vector used in the 1 + 1 + 2 
formalism. We can then construct a basis by introducing two spatial unit vectors in the 2-dimensional sheet (i.e. 
the screen), e° with / = 1,2 so that CjCja = Sjj and e^ria — ejUa = 0. It follows that h'^k^ = fc° + (1 + z)u°', 
hlk^ka = (1 + zf and /ige/a = e/b. 

The central equation governing gravitational lensing is the Sachs equation that derives from the geodesic deviation 
equation. Considering a geodesic in the bundle x"" = x'^+^°, where the vector can be decomposed as Co^^+X]/ f/^/- 
The geodesic deviation equation then takes the form 



dA 



where 72,/ = R'^i^cdk^k'^e-ia^^'^ ■ The linearity of the geodesic equation implies that it is related to the initial value of 
its derivative by a linear transformation ^/(A) = I?/(d^j/dA)o, so that the Jacobi matrix Vj relates the shape of the 
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cross-section of the bundle to the basis {ei,ej). For a bundle converging at the observer, ^/(O) = and the matrix 
Vj evolves according to the Jacobi equation 

^Vj = nfvi (B2) 

with I'/(0) = and (d2?//dA)o = 5/. Since the direction of observation is 6^ — (d^^/dA)o and the direction of 
the unlensed source 9^ — [Xg)/ Da{Xs), we conclude that the amplification matrix is -4/ (A) = {X) / D a{X) ■ This 
symmetric matrix is usually decomposed in terms of a shear (71,72) and a convergence k as 

Au^l^-""^^' , V (B3) 



2. The LTB case 



For an almost LTB universe, P = 0, (7° = 0, iTab = so that Mq = 0. It follows that the Riemann tensor can be 
decomposed as E^^d = pi?"^^ + Ei?"^^ + ni^^^d with 



P-K cd - 



p{uK[,h''^d]+h''[ch''d]) (B4) 



3 

Ei?°'cd = AuK[,E'^a]+4h^\cE'^d] (B5) 
Hi?°^d = 27j''''ui,Ha], + 2Tj,deU^-H''^'. (B6) 

The two first terms contain background and first order term while the third is purely first order. It follows that 

(1 + z)-^n-l = - [{inGp + £) 5j + 2E'ieiae''] - eiae'^m (v'^'^Hde + Vd'^H^) , (B7) 

where we have used the decomposition of k°-. Now, using the decomposition (|Aip for the electric and magnetic Weyl 
tensors (reminding that Nab — hab — naUb = e\e\ + e^e^), we obtain that 

(1 + zy^Uij = - [AnGp6ij + 2£a{e] + e}) + 2f,6e?e^j] + n^abe [nelje}^ - 2efje'j-^n'^ . (B8) 
At the background level, only the first term contributes so that 

7^/ = -47rGp(l + (B9) 

and the Sachs equation reduces to 

= -4^Gp(l + z)^Vj, (BIO) 

from which we deduce that Vj = f{X)lj. Since it is proportional to the identity matrix, it follows that the shear 
vanishes at the background level and we only have a convergence, exactly as in the Friedmann case. This was expected 
since for an observer at the center, the universe looks isotropic. The function / then satisfies 

^/ = -4^Gp(l + zff = -^Rabk^'k'f. 

The angular distance relates, by definition, the area of an object to the solid angle under which it is observed, 
dS^ — D\drt^ satisfies the same equation as / and has the same initial condition in [Da{0) = and dDA(0)/dA — 1] 
so that / — Da- This follows from the fact that if Dab is decomposed as 

^u^('\'l f-f'), (BID 

it can then be shown [U [H^l that the angular distance is related to the convergence by dlJ^i/dA = ODa- Now, 
derivating this equation and expressing the derivative of the convergence in terms of the shear on gets 50] that 

^Da = - Ul + ct| + ^Rabk''k''\ Da, 
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independently of the spacetime geometry. Now, in the particular case in which (T1/2 = 0, which is indeed the case 
at background level since 71/2 =0 we conclude that Dj^ satisfies the same equation as / and has the same initial 
condition in [Da{0) = and d£)A(0)/dA = 1] so that f = Da- 

Let us now turn to the perturbations. The expression (IBSP shows that the perturbation will enter both the shear 
and the convergence. Interestingly the shear vanishes at the background level. We can formally integrate the Sachs 
equation to get 

^"-i (^^2) 

Indeed the splitting TZij — TZ^'j'] + TZ\^] is ambiguous but can be used to obtain the expression of the shear since 
it vanishes at the background level. Note however that the gauge issue strikes only the convergence since the shear 
vanishes at background level. We conclude that, since 71 = {An — vA22)/2 and 72 = {A12 + ^2i)/2 (where the 
symmetrisation allows to get rid of the rotation that cannot be observed), we have 

^(7^ll+7^22) = '4nGp^2£aiel + el) 

^(7^ll-7^22) = -2£aiet ~ el) ^ V2£abPf 

i(7^l2 + 7^2l) - -2fa(e? + el) - V2£abPi'' - nSabeW^Pf, 



where 



We conclude that the shear is expressed only in terms of gauge invariant variables which is not the case of the 
convergence since it is non-vanishing at the background level. However K(A,n'') so that at a fixed A (redshift), one 
can extract the effect arising from the perturbation by taking a derivative in the sheet (since the convergence induced 
by the background is isotropic). 

We recover the standard FL expression when £ = at the background level. Here, the Sachs equation becomes 

^^vi = -\Hln,{i + z)^vi. 

It can be shown [s^, |5l| that the solution of this equation is exactly Da{,z) so that Vj = Da{z)Sj (but the general 
argument [13] ensures that it was expected). At perturbation level, Eab = DaDh^— ^A^hab = (9a9b$— ^A^Sab)/a'^ 
so that 

STlij = -{l + z)^[{47rGdp + d33'i>-A<P)di.j + 2DiDj<^>] 
= -{1 + z)^ [833^61 J + 2DiDj^], 

once the Poisson equation is used. 

The expression (|B12p gives the generalisation of the expression of the shear in terms of the perturbation variables 
of the 1+1+2 formalism for an observer seating at the center of a LTB universe. 

I 
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